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ABSTRACT 

The acceleration of relativistic jets from the Poynting to the matter dominated stage is 
considered. The are generally two coUimation regimes, which we call equilibrium and 
non-equilibrium, correspondingly. In the first regime, the jet is efficiently accelerated 
till the equipartition between the kinetic and electro-magnetic energy. We show that 
after the equilibrium jet ceases to be Poynting dominated, the ratio of the electro- 
magnetic to the kinetic energy decreases only logarithmically so that such jets become 
truly matter dominated only at extremely large distances. Non-equilibrium jets remain 
generally Poynting dominated till the logarithmically large distances. In the only case 
when a non-equilibrium jet is accelerated till the equipartition level, we found that the 
flow is not continued to the infinity but is focused towards the axis at a finite distance 
from the origin. 
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1 INTRODUCTION 



Collimated, Poynting dominated outflows are considered as 
a viable model for relativistic jets in active galactic nuclei 
(AGNs), microquasars and gamma-ray bursts (GRBs). For 
the relativistic outflows, the question of how the electro- 
magnetic energy is transformed into the plasma energy has 
no simple answer. In the non-relativistic case, the Poynt- 
ing flux is efficiently converted into the kinetic energy of 
the flow; an approximate equipartition is reached already 
at the Alfven point, where the toroidal magnetic fleld be- 
comes comparable with the poloidal field. Relativistic flows 
remain Poynting dominated even at the fast magnetosonic 
point. The reason is that in this case, the magnetic force is 
balanced not by inertia but by the electric force so that the 
plasma is only weakly accelerated by a small residual of the 
magnetic and electric forces. It turns out that in unconfined, 
nondissipative flows, the characteristic energy transforma- 
tion scale is inadmissibly large; such a flow is accelerated 

1/3 

only to the Lorentz factor of the order of 7max, where 7max 
is the Lorentz factor corresponding to the total transfor- 
mation of the electro-magnetic into the kinetic energy, af- 
ter which acceleration practically ceases, 7 ~ (7max In r)^''"^ 
(Tomimatsu 1994; Beskin, Kuznetsova & Raflkov 1998). 

The electro-magnetic energy could be more efficiently 
converted into the kinetic energy if the flow is collimated 
by an external medium. Such a configuration arises natu- 
raUy in gamma-ray bursts, where the relativistic jet from 
the collapsing stellar core pushes its way through the stellar 
envelope. In the accreting systems, the magnetically driven 



outflow from the rotating black hole could be collimated 
by the pressure of a slow (and generally magnetized) wind 
from the outer parts of the accretion disk. CoUimation and 
acceleration of externally confined, Poynting dominated jets 
has being studied extensively both numerically and analyti- 
cally, see the resent works by Komissarov et al. (2007, 2009), 
Narayan et al. (2008), Lyubarsky (2009, thereafter Paper I) 
and references therein. 

In the simplest case of the power-law external pressure 
distribution, p oc z"", the conditions for the flow accelera- 
tion and coUimation are the following (Paper I): 

(i) at AC > 2, the flow becomes asymptotically ra- 
dial and the acceleration is practically saturated at 7 ~ 
(7max/6^)^'''^, where O is the flnal coUimation angle, which 
itself is determined by the outer pressure distribution; 

(ii) at K ^ 2, the flow is accelerated until it ceases 
to be Poynting dominated; the shape of the flow line is 
paraboloidal, r ^ , where < 1 is determined by the 
outer pressure distribution. 

One sees that in the scope of ideal MHD, the electro- 
magnetic energy is efficiently converted into the kinetic en- 
ergy only if the flow is confined by the external medium 
with the pressure decreasing not too fast. However, one has 
to stress that at < 2, the flow Lorenz factor grows pro- 
portionally to the jet radius so that if the pressure decreases 
too slowly, so that the flow expands slowly, the acceleration 
rate would also be very low. In particular, if the surround- 
ing pressure goes to a constant, the flow becomes cylindrical 
and stops accelerating. 
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Till now only acceleration of the Poynting dominated 
flows has been addressed; the results obtained were just ex- 
trapolated to the energy equipartition stage. The transition 
from the Poynting dominated to the matter dominated stage 
has not been studied yet. An important point is that the flow 
could be considered as truly matter dominated if the ratio of 
the Poynting to the kinetic energy flux, a, becomes less than 
approximately 0.1. The reason is that only in this case, the 
shock jump conditions become close to those in the unmag- 
netized medium (Kennel & Coroniti 1984; Appl & Camen- 
zind 1988) so that the interaction of the jet with the sur- 
roundings occurs as in the non-magnetized case. At a > 0.1, 
only weak shocks are possible therefore the flow pattern, 
which arises when such a jet is decelerated in the ambient 
medium, significantly differs from that for the purely hy- 
drodynamic jet (Komissarov 1999). Simulations show that 
already at cr « 0.01, the flow pattern differs significantly 
from the purely hydrodynamic one (Leismann et al 2005). 
In the GRB context, Mimica, Giannios and Aloy (2009) and 
Mimica and Aloy (2009) show that even a moderate mag- 
netization of the ejecta could have a profound effect on the 
properties of the internal shocks as well as on dynamics of 
the deceleration thus affecting both the prompt and the af- 
terglow emission. 

Here we study the transition of the flow through the 
cr ~ 1 domain. We address only the case k ^ 2 when such a 
transition could occur at all. The paper is organized as fol- 
lows. In Section 2, we present the asymptotic equations de- 
scribing relativistic, magnetized flows at large distance from 
the origin. In Section 3, we find solutions to these equations 
in the case k < 2. The case k = 2 is addressed in Section 3. 
Conclusions are presented in Section 5. 



2 THE JET STRUCTURE IN THE FAR ZONE 

In Paper I, we presented asymptotic equations describing 
the relativistic, magnetized flow at distances much larger 
than the light cylinder radius. Now we shortly outline the 
relevant results. As usual, the magnetic field is conveniently 
decomposed into the poloidal and toroidal components, B — 
Bp + B^(f), the poloidal field being expressed via the flux 
function 



Bp = -V* X 

r 



(1) 



We use cylindrical (r, 0, z) coordinates; the hat denotes unit 
vectors. The distribution of the mass flux at the inlet of 
the flow is described by the function rji}^) defined by the 
continuity equation 



47r/9Wp7 

Br, ' 



(2) 



where p is the plasma density in the lab frame, 7 the flow 
Lorentz factor, Vp the poloidal velocity. The conserved total 
energy flux is presented as 



7- 



■n 



(3) 



where is the angular velocity of the flied line. In this 

expression, the first term is the kinetic energy whereas the 
second one is the Poynting flux. 



The energy integral p is determined from the condition 
of smooth passage of the flow through the singular points. In 
the far zone, it should be considered as a given function. An 
important point is that close to the axis, the energy integral 
has the universal form 



M(*) = Tin ( 1 + ^ 



JinV. 

2^2 ' 



(4) 



where 7in is the Lorentz factor at the inlet of the flow. In 
this expression, the first and the second term describe the 
kinetic and the Poynting energy flux, correspondingly, at the 
inlet of the flow. Note that the Poynting flux goes to zero 
at the axis therefore the flow is Poynting dominated only at 
* > *. 

The structure of the flow is described by the transfield 
and Bernoulli equations. As an unknown function, one can 
conveniently use the shape of the magnetic flux surface, 
r{z,'ii). For coUimated flows, r <C z, the transfield force 
balance equation in the far zone could be written as 



+ 



1 



1 - 



+ 



27i 

^l 7 



(5) 



The Bernoulli equation is reduced, beyond the fast magne- 
tosonic point, to 



dr 
9* 



(6) 



Ifere we are interested in outflows subtending a finite 
magnetic flux therefore Eqs. ((5]) and ((6)1 should be solved 
at ^ *]/ ^ ^0- We assume that 'I'o ^' so that the main 
body of the flow is initially Poynting dominated. If the flow 
is confined by the pressure of the external medium, pcxt(2:), 
the pressure balance condition should be satisfied at the 
boundary: 



^(m ^ 7) 



= 87rpext(2). 



(7) 



The boundary condition at the axis is r{'^ = 0) = 0. 

There are generally two different regimes of coUimation. 
At the condition 



dz2 



1 



(8) 



one can neglect the term with the derivative in z and write 
the transfield equation (O as an ordinary differential equa- 
tion 



M 1 + 



7^ 



- 27in 



n^r^ + 'yl, , d r7(/i-7) 



-(m-7) 



9* ^7 



(9) 



This equation describes in fact cylindrical equilibrium, in 
which case the residual of the magnetic hoop stress and 
the electric force is counterbalanced by the pressure of the 
poloidal field. The corresponding coUimation regime is called 
equilibrium because in this case, the structure of the jet at 
any distance from the origin is the same as the structure of 
an appropriate equilibrium cylindrical configuration. For a 
smoothly expanding jet, the condition ((Sjl is reduced to 



< z; 



(10) 
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Note that neglecting the second derivative in the transfield 
equation ((5}, one looses solutions. These lost solutions just 
describe oscillations of the flow with respect to the equilib- 
rium state satisfying Eq. ((Ojl . Therefore if the condition (|10p 
is fulfilled, one can anyway use Eq. (|9} in order to find the 
overall expansion of the jet. 

The transverse equilibrium implies that in the proper 
plasma frame, the toroidal and poloidal magnetic fields are 
comparable. Bp ~ B'^. Transforming to the lab frame, one 
gets 7 ~ B^/Bp. Taking into account that the toroidal field 
is wound up from the poloidal one so that B^ « QrBp, one 
concludes that in the equilibrium flow, 7 ~ Qr. This esti- 
mate is confirmed by explicit solutions (Tchekhovskoy et al 
2008; Komissarov et al 2009; Paper I; Beskin & Nokhrina 
2009). Now one sees that the condition (|10p in fact implies 
that the flow is in causal connection, i.e. the proper propa- 
gation time 2/7 exceeds the time r necessary for a signal to 
cross the flow so that the flow has enough time in order to 
settle into transverse equilibrium. 

If the condition opposite to (|10[) is fulfilled. 



(11) 



which anyway could happen only far enough from the axis, 
the term with the second derivative becomes dominant and 
the transfield equation is reduced to 



2anr—-—^ 



7f 



(12) 



This equation could be directly obtained neglecting 
the poloidal magnetic field and the azimuthal velocity 
(Lyubarsky & Eichler 2001). Such a fiow may be con- 
ceived as composed from coaxial magnetic loops. We call 
the corresponding collimation regime non-equilibrium. Non- 
equilibrium flows are causally disconnected or marginally 
connected therefore the magnetic loops generally do not 
shrink even though the poloidal held pressure is negligibly 
small (see, however, sect 4). 

In Paper I, the structure of the jet was found under the 
condition that the main body of the flow remains Poynt- 
ing dominated. Here we relax this condition and study the 
transition to the matter dominated flow. We consider the 
jet with a constant angular velocity, Q, = const and homo- 
geneous injection, 77 = const; 7in = const. In this case, one 
can conveniently use the dimensionless variables 



X ^Q.r; Z = Q,z. 



(13) 



We also assume that the energy integral, ^^{^), is described 
by a linear function @ throughout the flow. The last is a 
good approximation for jets with a constant angular velocity 
(Komissarov et al 2007, 2009; Tchekhovskoy et al 2008). We 
assume that the external pressure decreases as a power law 



P = 



Po 



Po 



Z« ' 67r 
The normalization coefficient /3 is chosen as in Paper I 



(14) 



3 EQUILIBRIUM JETS 

The Poynting dominated jet is collimated in the equilibrium 
regime if «: < 2. Then the jet expands as X (x whereas 
the Lorentz factor increases as 7 ~ X (Tchekhovskoy et al 



2008; Komissarov et al 2009; Paper I; Beskin & Nokhrina 
2009). Now we address the transition of the equilibrium jet 
from the Poynting dominated to the kinetic energy domi- 
nated flow. 

In the equilibrium regime, the flow is described by a pair 
of ordinary differential equations ^ and Q for the trans- 
verse structure of the jet, the dependence on Z entering only 
via the boundary condition ([7]). Introducing the variables 



X 7 * 

7in 7in * 

one reduces Eqs. (O and © to the dimensionless form 



= 1 - 



2{i + e){s~r)' 



r)- 



(15) 

(16) 
(17) 



r ds 

ds 2(s 

Solutions to these equations were analyzed in detail in Paper 
I. Any solution describes the transverse structure of the jet 
at some distance from the origin. Near the axis, ^ ^ 1, the 
solution is 



i+ce-, 



(18) 



where C is a constant. If C > 0.38, the solution goes to a 
Poynting dominated flow far from the axis, ^ ^ 1. In the op- 
posite case, C < 0.38, the solution describes the transverse 
structure of the matter dominated jet. 

For any constant C, the solution to Eqs. (I16|l and 
(|17|) could be easily found numerically by the Runge-Kutta 
method. In order to flnd the full structure of the jet, one 
has to find an appropriate constant C for any distance from 
the origin, Z. We are looking for a solution satisfying the 
outer boundary condition ((7|, which is written in the new 
variables as 



where 

7m ax 



s = l-t-*o/* 



P 7in7. 



(19) 



(20) 



is the maximal achievable Lorentz factor of the fiow, which 
is just Michel's magnetization parameter (Michel 1969) . At 
any Z, one has to find a constant C such that the solution 
to Eqs. (I16|) and (|17[) satisfies both the condition (|18p and 
the condition (|19[) . This could be easily done by bisection. 

As an example, we presented in Fig. 1 the structure of 
the jet confined by the pressure p oc z^^"^ . One sees that ini- 
tially the flow expands according to the Poynting dominated 
scaling X oc Z'^^'^ — Z^^^. A cylindrical, moderately mag- 
netized, cr ~ 1, core is formed within the jet at this stage, 
as it was discussed in Paper I. When the bulk of the flow 
ceases to be Poynting dominated, the jet begins to expand 
faster than at the Poynting dominated stage. Note that in 
the logarithmic plot, it looks as if the jet inflates however, 
the jet collimation angle, X/Z, still decreases but slower 
than at the Poynting dominated stage. The weaker collima- 
tion of the flow in the matter dominated stage has already 
been noticed by Komissarov et al (2009). In Fig. 2, we show 
evolution of the Lorentz factor and of the magnetization 
parameter along selected flow lines. Note that according to 
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the general equilibrium scaling 7 ~ X, the Lorentz factor is 
larger at the periphery of the flow. One sees also that even 
though a continuously decreases with the distance, the a 
decreasing rate drops down after the flow reaches equipar- 
tition therefore the flow becomes truly matter dominated 
(cr < 0.1) only at extremely large distances. 

Let us now find analytically asymptotic solutions de- 
scribing equilibrium flows in the Poynting and matter domi- 
nated domains, correspondingly. With this purpose, one can 
conveniently use cr as a variable instead of s: 



Then Eqs. (fT6|) and (fTfll are reduced to 
da 



(21) 



(22) 



r 



= {1 + if - r^)cr - (r - 1)^(23) 



The condition at the axis is cr(0) = 0. The outer boundary 
of the jet, is defined by the condition 

s(Co) = 1 + *0/* ~ *o/* = 7max/7m; (24) 

then Eq. ([21]) yields 

-(Co) = - 1. (25) 

The boundary condition (I19|l is written in the new variables 



Co 



(26) 



Let us first find the solution for the Poynting dominated 
part of the jet, cr ^ 1. Note that since the Poynting flux goes 
to zero at the axis of the flow (see Eq. Q), the condition 
cr 3> 1 could be met only at ^ 1. Moreover, we can take 
r 3> 1 in this range because the Poynting dominated flow is 
accelerated so that eventually the Lorentz factor of the flow 
exceeds the initial one. Then Eqs. ([55} and ((25)) are reduced 
to 



C3 rfcr 2 

di 



The solution to the second equation is 



(27) 
(28) 

, (29) 

^l+cxp(-Z)C2) 

where D is a constant. Since Eqs. (|27|l and (|28p are valid 
only far from the axis, one has to solve them in the limit 
^ ^ 1; this yields 



r = 



r = C; <y = Ai- 



(30) 



where yl is a constant, which is found from the boundary 
condition (|26[) as 



3 



(31) 



Note that the ^ 3> 1 solution is independent of the constant 
D therefore the structure of the flow far from the axis is 
uniquely determined by the outer boundary condition (|26|) . 



The condition at the axis of the flow, Eq. (|18|) . does not place 
any restriction on the structure of the Poynting dominated 
flow at ^ > 1. 

The expressions (I30p and pi|l describe the internal 
structure of the equilibrium jet not too close to the axis. 
One sees that the Lorentz factor at any point of the flow is 
proportional to the cylindrical radius of the point, 



7 ■ 



(32) 



so that in any cross-section of the jet the Lorentz factor 
increases outwards whereas at any flux surface it increases 
with the distance so far as the flow expands. This is the 
general property of equilibrium Poynting dominated jets 
(Tschekovskoy et al 2008, Paper I, Beskin & Nokhrina 2009). 
The jet radius may be found by substituting Eqs. (|30|l to 
the condition p5|) and taking into account that a ^ 1; this 
yields 

Xo = 7ineo = ' . (33) 

This expression was already obtained, by different methods, 
in (Tschekovskoy et al 2008, Komissarov et al 2009, Paper 
I)- 

The presented solution remains valid while o remains 
large, i.e. outside the boundary defined from Eqs. (|30|l and 

m 



(34) 



Near the axis, a cr « 1 core is formed (Paper I). Comparing 
Eq. (|33p with Eq. (I34p . one sees that the fraction of the jet 
volume occupied by the core grows with the distance from 
the origin therefore eventually the whole jet ceases to be 
Poynting dominated. This happens at the distance 

^cquipait = I 5 I ■ (35) 



The full solution presented in Figs. 1 and 2 confirms this 
scaling. 

Let us now solve Eqs. (|22p and (|23|l in the zone Z ^ 
^oquipart, whcrc the Poynting fiux is already small as com- 
pared with the kinetic energy fiux, cr ^ 1. Inspecting Eq. 
(|22p . one sees that in order for cr to be small, P should be 
close to unity at ^ ~ 1. Moreover, comparing the left-hand 
side of Eq. (|23p with the first term in the right-hand side, 
one sees that there should be P — 1 ~ a. Then one can ne- 
glect (r — 1)^ with respect to a and reduce Eqs. (|22p and 

to 



2a; 



= 2crr 



r d£_ 
The solution to Eq. ((Sejl is 



(36) 



(37) 



(38) 



1 + e' 
where C is a constant. 

One sees that a goes to a constant, cr = C, at a large ^, 
which agrees with the general analysis presented in Paper I. 
As this solution is obtained under the condition ct ^ 1, one 
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lg(^/7in) 

Figure 1. The flux surfaces in tlie equilibrium jet; k = 1.5; 7max/7in = 200. The equipartition radius is defined by Eq. 11351 . Dotted line 
shows the Poynting dominated asymptotics 11331 whereas dashed line shows the matter dominated asymptotics II49I I. The discrepancy 
between the last and the numerical solution is attributed to the fact that even at the distances as large as lO^Zoquiparti o" is still not 
very small (see Fig. 2). 



2.5 




-^equipart) 

Figure 2. Evolution of the Lorentz factor (dashed) and of the ratio of the Poynting to the kinetic energy fiux (solid) along the boundary of 
the jet (thick lines) and along the flow line marked by asterisks in Fig. 1 (thin lines). One sees that, according to the general equilibrium 
scaling 7 ~ Hr, the Lorentz factor is larger at the periphery than within the jet. One also sees that at Z 3> ^oquipart, o" becomes 
independent of r and decreases with z very slowly; this is confirmed by the asymptotic result of Eq. JSOj. 
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concludes that there should he C <^ 1. With account of Eq. 
Eq. dSZll yields 



r = 1 + c Hi + C) 



1+e 



(39) 



Expanding Eqs. p8|) and (|39|l in small ^, one sees that the 
constant C in this solution is the same that in Eq. (|18p . 
Recall that the solution (I38II and (1391) was obtained under 



the assumption (F — 1) ^ a ^ 1, which implies C <^ 1 



and ln(l <C 1/VC The last condition shows that the 

solution (|38p and (|39|) is valid from ^ = to a large ^ but 
not too large, 



5 <C exp 



2VC 



(40) 



Let us now find the solution for an unrestrictedly large ^ 
which could be matched with the solution (|38p and (|39|l . 

One can easily find the solution to Eqs. (|36[) and (|37| 
at the condition 



Then Eq. ((22l) is reduced to 



< cr; 



(41) 



(42) 



which yields a = const. This solution is matched with the 
solution (|38p if o" = C so that the solution (I38p could be 
continued to an arbitrary large ^. In the same limit, Eq. 
(1231) is reduced to 



r 



= 2(7 = 2C; 



(43) 



which yields F = Ci^^'^, where Ci is a constant. In the 
region (I40p . this solution is reduced to F = Ci(l + 2Cln^), 
which is smoothly matched with the solution (|38p provided 
Ci = 1 — C. Now one can write 



(44) 



Recall that we solved the equations under the condition (|4ip . 
The function (|44p . together with the function ct = C ^ 1, 
satisfy this condition at ^ ^ l/VTJ. The Lorentz factor of 
the flow is given by the expression (|39} at In ^ ^ 1/ x/C* and 
by the expression (gj)) at ^ >• 1/y/C, these two expressions 
being smoothly matched in the region 



c 



(45) 



One finally concludes that Eqs. (I38p . (|39p and (|44p represent 
the full solution for the transverse structure of a low-cr jet. 

The dependence of the jet structure on the distance 
from the origin, Z, enters via the constant C, which is found 
from the boundary conditions (|25p and (|26p . They could be 
written, substituting C for (t(5o) and Eq. (|44p for F(5o), as 

^max 



7in 




■c/2 



Taking into account that C <^ 1, one gets 

ln(7max/7in) 



C - 



21n^o 



(46) 



(47) 



(48) 



27iiigo In go 

7maxln(7max/7in) 



2/k 



(49) 



The last two equations provide dependence of the constant 
C on Z thus closing the solution for the low-cr part of the 
jet. 

Eq. (|49p describes the shape of the flow lines at the 
a <^ 1 stage, see Fig. 1. Taking into account that this so- 
lution is valid at k < 2, one sees that the collimation angle 
continuously decreases, d^o/dZ — > 0, so that the jet becomes 
asymptotically cylindrical. According to Eq. (|38p . C is equal 
to CT at g ^ 1, i.e. not too close to the axis. Combining Eqs. 
and (|49p . one finds an estimate 



'Tie » 1) 



ln(7max/7in) 



(50) 



2 ln(7max/7iii) + (f^/Z) ln(Z/Zcquipart) ' 

Note that a is constant across the jet at g ^ 1; the full 
solution presented in Fig. 2 confirms this asymptotic result. 
According to Eq. (|50p . a becomes of the order of unity at 
Z ~ ^cquipart, SO that the obtained asymptotics is roughly 
matched with the asymptotics (|30l33p for the Foynting dom- 
inated jet. An important point is that at Z > Zcquipart, o" 
decreases however extremely slowly. For example, one sees 
in Fig. 2 that for the chosen parameters of the jet, a de- 
creases only to 1/3 at the distance as large as 10*.Zequipart. 
The fiow becomes truly matter dominated, a < 0.1, only at 

Z > 

■^cquipart (7n.ax/7in)*/'' . (51) 

The results of this section are valid only for equilib- 
rium jets; the corresponding condition is given by Eq. pop . 
It follows from Eq. (|33p that the Foynting dominated jet 
is collimated in the equilibrium mode if the outer pressure 
decreases slowly enough, «: < 2. At the matter dominated 
stage, the fiow is collimated slower. One sees from Eq. (|49p 
that only if k < 1, the flow at this stage satisfles the con- 
dition (|10p till the infinity. At 1 < k < 2, the fiow expands 
faster than Z oc therefore eventually the condition for 
the equilibrium collimation is violated. This occurs at the 
distance 



2— 



(!) 



Tmax 



(52) 



The transition from the equilibrium to the non-equilibrium 
regime could be studied only numerically. 



4 NON-EQUILIBRIUM JETS 

In this section, we study transition through cr ~ 1 in non- 
equilibrium jets. The jet is efficiently accelerated in the non- 
equilibrium regime only ifK = 2;/3<l/4 (Faper I). At 
K < 2 the jet is accelerated in the equilibrium regime. In 
the case n — 2; P > 1/4, the acceleration occurs in the 
intermediate regime when the jet is not in the transverse 
equilibrium however one cannot neglect the poloidal field; 
practically this intermediate regime is close to the equilib- 
rium one. At K > 2 the jet is in the non-equilibrium regime 
however, the acceleration is saturated at a terminal Lorentz 
factor, which is generally less than 7max. 

Of course if the initial Foynting fiux is not too large, 
the fiow could be accelerated till the ~ 7max even if «; > 2. 
For example, even though the non-confined fiow is acceler- 
ated only till the terminal Lorentz factor 7 ~ (7max Inr)^^^ 
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(Tomimatsu 1994; Beskin et al. 1998), the equipartition is 
reached close enough to the axis, where the Poynting flux 
does not exceed this limiting value (Lyubarsky & Eichler 
2001; Tchekhovskoy, McKinney & Narayan 2009a). Another 
example is the flow confined by the external pressure with 
K something larger than 2. Then the flow is accelerated to 
the terminal Lorentz factor 



7t ~ (7max/0 ) 



2x1/3. 



(53) 



where O is the flnal coUimation angle, which is determined 
by the outer pressure distribution (Paper I). Such a flow 
could reach equipartition only if 7max ^ 7t. An important 
point is that this is possible only if 7max6 ^ 1 so that such 
an efficient transformation of the Poynting to the kinetic en- 
ergy anyway occurs only in causally connected, equilibrium 
flows. One should note that even if 7max > 7t, so that the 
bulk of the flow stops accelerating still being Poynting dom- 
inated, the flow in the boundary layer is accelerated till the 
equipartition provided the confining pressure decreases; the 
width of this boundary layer is determined by the condition 
of the causal connection with the boundary (Paper I, see also 
fig. 3a in Tchekhovskoy, Narayan & McKinney 2009b). We 
do not address all these cases in this paper and concentrate 
on the conditions permitting unrestricted acceleration up to 
7 ~ 7max. The non-equilibrium jet could be accelerated till 
7 ~ 7max for any 7inax only iffc = 2,/?<l/4. Here we study 
this case. 

The non-equilibrium jet is described by the transfield 
equation in the form (|12|) , which should be supplemented by 
the Bernoulli equation Eliminating 7 from these equa- 
tions and making use of Eq. Q for ^, one gets a single 
equation, which looks in the dimensionless coordinates X 
and Z as 



d 



X 



ax 

df-i 



X dfi 2^i^ 



df-l 



X 



(54) 



This equation is invariant with respect to the transformation 
fj. sjj., X — > s™X, Z s^'^^Z, where s and m are 
arbitrary numbers, therefore one can look for a self-similar 
solution in the form 



X = M'"f/(C); C = 



(55) 



With this ansatz, Eq. (|54|l is reduced to an ordinary differ- 
ential equation for U : 



[(2m - 1)U - 2(1 + m)CU'] ^ U" 
= 4(1 + mfCpU' + CUU" - CU% 



(56) 



where the prime denotes derivative with respect to (. 

The general solution to Eq. (|56p is presented in Ap- 
pendix. We are interested in a solution satisfying the bound- 
ary condition (O with the pressure distribution p4p . k — 2. 
With the ansatz (|55p . this boundary condition is written, 
with the aid of the Bernoulli equation ((6}, as 



(2m - l)f/ - 2(m -M)C(7' 



(57) 



Here we used Eq. @ at = 'I'o 2> 'P. So we are looking 
for a function satisfying both Eq. (|56p and Eq. (|57p . The 
solution to Eq. (|57|) is written as 



U = 



where ("0 is a constant, 
2m- 1 

k = 



2(1 



(58) 



(59) 



One can check straightforwardly that this function satisfies 
also Eq. (|56p provided 



(60) 



fc = ^ (1 + - 4/3) 



Note that 1/2 < fc < 1. This solution exists only if /? < 1/4, 
which agrees with the conclusion in Paper I that at k = 2, 
the flow is in the non-equilibrium regime only at /3 < 1/4. 

Returning to the physical variables (|55|) . one finds the 
shape of the flux surfaces (see Fig. 3) 



X = 



3/2 



(61) 



The Bernoulli equation (O provides the expression for the 
Lorentz factor of the flow 



<,i-* + 2c>-'" 

At ( <^ (0, one comes to the scaling 

XocZ''; ■ycxZ^''' (63) 

obtained earlier (Komissarov et al 2009; Paper I) for a 
Poynting dominated jet at k = 2, /3 < 1/4. The expressions 
(|61|) and H62|) generalize this solution beyond the Poynting 
dominated domain. Note that when /3 goes to 1/4 from be- 
low, k goes to 1/2 so that one comes to parabolic flow lines, 
X oc v'^, in which case one has to take into account the 
neglected poloidal field. In the opposite limit, /3 = 0, one 
comes to an unconfined radial flow, k = 1. 

The constant ^0 should be found by matching this so- 
lution with the near zone solution. If the flow is not coUi- 
mated at the light cylinder, X ~ 1, Z ~ 1, there should be 
Co 

be Poynting dominated at the distance 



7max'^'^ • Then one concludes that the jet ceases to 



7 A 3/[2(l-fe)l 

'^cquipart ~ S07max 



, 1/(1-'=). 

/max 1 



(64) 

which agrees with the estimate presented in Paper I. 

One sees that when approaches (0, the flow converges 
to the axis and the Poynting flux is converted into the kinetic 
energy, 7 ^. Note that at the Poynting dominated stage, 
C ^ Co, the flow was causally disconnected, 



X7 
~Z' 



= W^>1. 



(65) 



When it approaches equipartition, the bulk Lorentz factor 
does not grow any more, the causal connection is estab- 
lished and the magnetic loops squeeze the flow. Recall that 
all the results of the present section are obtained at the as- 
sumption that one can neglect the poloidal field; the corre- 
sponding condition is given by Eq. I|lip . If the flow remains 
axisymmetrical, it stops converging in the region X y/Z 
where the pressure of the poloidal field becomes significant. 
However, one can expect that due to the kink instability, 
different magnetic loops could come apart in the converging 
flow forming an irregular field structure, which could trigger 
dissipation processes. 
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Figure 3. The flux surfaces in the non-equilibrium jet; ft = 2; /3 = 0.2; k = 0.72; 7max = 100. Dotted line shows the Poynting dominated 
asymptotics X = = Z'^-'^^. Dashed line shows the boundary X = yfZ within which the poloidal field could not be neglected therefore 
the presented solution is valid only to the right of the dashed line. 



5 CONCLUSIONS 

In this paper, we studied the acceleration of Poynting domi- 
nated jets to the matter dominated stage. Even though effi- 
cient transformation of the electromagnetic into the kinetic 
energy is possible in principle in the scope of ideal MHD, 
the conditions for such a transition are not trivial. Namely, 
the flow should be confined by external pressure, which de- 
creases with the distance but not too fast. It was shown 
in Paper I that the flow is accelerated until it ceases to be 
Poynting dominated only if k ^ 2 (for the power law con- 
fining pressure distribution 113}). In the opposite case, the 
flow practically stops accelerating after it reaches a terminal 
Lorentz factor of Eq. (|53[). which is generally less than 7max. 
If 7t > 7max, the flow behaves essentially as in the k < 2 
case. 

Especially restrictive is the fact that the acceleration 
zone spans a large range of scales so that one has to ensure 
that the conditions for acceleration are fulflUed along all the 
way. If the confining pressure begins to decrease faster than 
oc z~'^ or drops down abruptly, which happens for example 
when the GRB jet escapes from the progenitor star, the 
acceleration and collimation cease and the jet propagates 
further out preserving the acquired collimation angle and the 
Lorentz factor. On the other hand, if the confining pressure 
stops decreasing, the flow becomes cylindrical however the 
acceleration is terminated. 



The flow is efficiently accelerated in the equilibrium 
regime, i.e. if at any distance from the origin, the structure 
of the flow is settled into the structure of an appropriate 
cylindrical configuration. At the Poynting dominated stage, 
the equilibrium flow is accelerated as 7 ~ S7r (Tchekhovskoy 
et al 2008; Paper I, Beskin & Nokhrina 2009) therefore the 
faster the flow expands, the faster it is accelerated. On the 
other hand, the condition (|10p implies that the flow is in the 
equilibrium only if it remains within the parabola X ~ ^fZ 
(which in fact ensures that the causal connection is main- 
tained accross the flow) therefore the acceleration rate is 
maximal for the parabolic flow. The flow expansion is de- 
termined by the distribution of the confining pressure. The 
faster the outer pressure decreases, the faster the jet ex- 
pands. Taking into account Eqs. (|10p and (|33|l. one sees that 
at «: < 2 the flow is in the equilibrium regimqj. 

The fastest acceleration regime is achieved when k, goes 
to 2. At K = 2, the flow is in the equilibrium regime only at 
/3 S> 1. The case /3 ~ 1 is an intermediate between the equi- 
librium and non-equilibrium regimes. At the Poynting domi- 
nated stage, the properties of the k = 2 flow are very similar 

^ If K is a bit larger than 2 and /3 > 1, the equilibrium conditions 
still could be met not too far from the source. In this limited 
region, the flow behaves as an equilibrium one, i.e. it expands 
and accelerates according to the equilibrium scalings (1321 and 
II33I I. see paper I. 
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to those of the equihbrium flow at /9 > 1/4 (Komissarov et 
al. 2009, Paper I). In this case, the jet has a paraboUc form 

(66) 

whereas the Lorentz factor grows as 7 = X. The energy 
equipartition is achieved in such a flow at the distance 
■^equipart ~ 7max- For a slower decreasing external pressure, 
K < 2, the flow remains Poynting dominated at even larger 
distances, see Eq. (I35|) . 

GRB jets are known to have Lorentz factors of at least 
a few hundreds (Lithwick & Sari 2001). In the fastest accel- 
eration regime, one could reach these Lorentz factors if the 
size of the acceleration region is ~ 10^(7/300)^ of the light 
cylinder radii. Taking into account that the characteristic 
light cylinder radius of a rapidly rotating black hole of a 
few solar masses is ~ 10^ -^ 10*^ cm, one sees that the above 
estimate of the acceleration zone is compatible with the size 
of the progenitor star. The problem is that the observations 
of the afterglow light curves as well as the burst statistics 
evidence for 76I ~ 10 30 (e.g. Tchekhovskoy et al, 2009b) 
whereas an efficient transformation of the Poynting to the 
kinetic energy occurs only if 7S ~ 1. Highly coUimated but 
causally disconnected jets could be formed if the confining 
pressure decreases something faster than 2"^ (Paper I) how- 
ever, such jets remain Poynting dominated so that magnetic 
dissipation is necessary in order to utilize the energy of the 
outflow. 

Observations of jets in AGNs and microquasars evi- 
dence for Lorentz factors from a few to a few dozens (e.g. 
Cohen et al 2003; Mirabel & Rodriguez 1999). In order 
to achieve these Lorentz factors in the fastest acceleration 
regime, one needs the size of the conflnement zone of only 
.^equipart ~ 100 lOOOrg provided the black hole is rapidly 
rotating. The wind from the accretion disk could serve as 
the confining medium up to distances of the order of the ex- 
ternal disk radius, which could be that large. It is not clear 
whether the wind from the disk could provide the confin- 
ing pressure decreasing not faster than r~^. If the pressure 
decreases faster, the jet remains Poynting dominated. An 
important point is that even if the necessary conditions are 
fulfilled, the flow could reach only an equipartition state 
but not a true matter dominated stage, cr < 0.1. In this pa- 
per we have shown that when the equilibrium jet ceases to 
be Poynting dominated, the collimation angle still decreases 
even though slower than at the Poynting dominated stage. 
However, a decreases only logarithmically so that ct ~ 0.1 is 
achieved only if the confining medium is extended beyond 

the distance .^equipart7max ~ 

10«^ 10^ (7,nax/ 10)^9 (see Eq. 
(|5ip . we assumed that 7i„ ~ 1), which seems to be inappro- 
priately large. The fact that without magnetic dissipation, 
jets could not become true matter dominated, has impor- 
tant implications for the interaction of the ejected material 
with the surroundings (Leismann et al 2005; Mimica & Aloy 
2009; Mimica et al. 2009). 

Since the causal contact is maintained across equilib- 
rium jets, one has to worry about the kink instability, which 
could significantly disturb or even destroy the regular flow 
structure. However, last studies reveal (Tschechovskoy et 
al. 2008) that in Poynting dominated outflows, the poloidal 
field is very close to uniform (and is exactly uniform for the 



chosen here simple expression Q for ^i(^) and — const, 
rj = const); in this case the growth rate of the kink instabil- 
ity goes to zero (Istomin & Pariev 1996, Lyubarskii 1999). 
It is possible that the instability could develop in spite of 
the low growth rate because the jet acceleration zone is very 
large but in order to clarify the question, more careful inves- 
tigation of the transverse structure of the jet is necessary. 
Note also that even if the instability turns out to be sup- 
pressed in the Poynting dominated stage, it could develop 
at the moderately magnetized stage when the poloidal field 
is concentrated to the axis of the flow (Paper I, Beskin & 
Nokhrina 2009). In any case the impact of the instability on 
the jet structure should be studied only with 3D numerical 
simulations. 

The acceleration in the non-equilibrium regime is gen- 
erally not very efficient so that the fiow Lorentz factor could 
not significantly exceed a terminal value determined by the 
parameters of the fiow and of the surrounding medium (Pa- 
per I). Only in the specific case k = 2, /3 < 1/4, the jet 
is accelerated in non-equilibrium regime till the equiparti- 
tion level. We showed that after such a fiow reaches rough 
equipartition, it sharply converges to the axis and the en- 
ergy is efficiently transferred to the plasma. It is not clear 
what happens to this "collapsing" fiow; one can expect that 
a sort of a "hot spot" appear in such a flow, which resem- 
bles that formed by hydrodynamical recoUimation of a rel- 
ativistic outflow (Levinson & Bromberg 2008; Bromberg & 
Levinson 2009). 

In any case, we have not found, in the scope of ideal 
MHD, a possibility for a smooth acceleration of a Poynting 
dominated flow to the matter dominated stage at a reason- 
able scale. In equilibrium jets, a decreases too slowly after 
the flow ceases to be Poynting dominated. In the only case 
of efficiently accelerated non-equilibrium jet, the flow is not 
continued till the inflnity but in fact collapses. 

In this paper, we have addressed the transition to the 
matter dominated stage both in the equilibrium and in the 
non-equilibrium flows. The intermediate case, k = 2; /3 ~ 1, 
is beyond the scope of our analytic approach. At the Poynt- 
ing dominated stage, the behavior of the intermediate flow 
is similar to the behavior of the equilibrium flows (Komis- 
sarov et al 2009; Paper I). In order to find what happens to 
such a fiow at the a < 1 stage, one has to solve Eqs. ((5)1 and 
([B} numerically. 
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Substituting the expressions HA4p and ^A5\i into Eq. (|55|l . 
one gets a set of non-equilibrium jet configurations depend- 
ing on three parameters, m (or k), c and ^o- The Bernoulli 
equation ([6]) provides the expression for the Lorentz factor 
of the flow: 



7 = 



1 + fc 



(A7) 



(2/3)(l-fc)- 

One sees that there should be <j!> > 1 — fc in order for 7 to be 
positive. Note also that 7 ^ /i (i.e. (t — > 0) only if <^ — > 00. 
The solution (|58|l , (|62|l is obtained from the general solution 
(fUl) . (fXS)) and d^ll) at c = 0. 



APPENDIX A: 

The general solution to Eq. (|56|l is found taking into ac- 
count that this equation is invariant under the transforma- 
tion C, — > sC,, U — > sU . Therefore we can reduce the order of 
the equation using the substitution 

/ = [//C; = (Al) 
/ dC, 

Then one gets a linear first order equation. 



m) 



\2 



^1 ^^^-^^ .2 , 4(l + .,„, 
2^^-'^^^=^ +[3_2(l + m)0]3' 

which is immediately solved to yield 

^^(^)-T^((^-(r4^)^ (^3) 

where c is a constant. Substituting then /(</>) to the sec- 
ond equation (|A1|) . one gets a diff'erential equation for C,{<j>), 
which is easily solved. 

Finally one gets the solution to Eq. (|56p in the para- 
metric form: 

^w = Coy ^ 3), V (i-'=)[i-(i-'=)-c^i (A4) 

l-(l-fc)^e^ 

(((?!>) = Co|l - 010 (i-'=)[i-(i-'=)"c"l (A5) 
where C,o is a constant, 

fci.2 = lT-^. (A6) 
1 ± c 



